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ABSTRACT 


The polarization rotation is studied for radio signals received 
from an artificial satellite in the ionosphere. A first-order analysis 
is used which leads to values for the integrated electron density in 
the ionosphere. These values are compared with those calculated by 
other means and with those obtained by other investigators. A new 
parameter, the rate of polarization null position, is introduced as a 
method for obtaining the electron density at satellite altitudes. The 
methods employed in the analysis axe outlined and additional means 


of improving the analysis are discussed. 
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LIST OF SYMBOLS 


(All units rationalized MKS where not specified) 


electron content in a column of height hg and ! square meter 


in area 
frequency in cycles per second 
height of the satellite in meters 


velocity in meters per second 


electric field of an electromagnetic wave 
electronic charge 

magnetic field density in webers. per meter” 
permeability of free space 

magnetic field intensity in amperes per meter 
force in Newtons 

mass of an electron 

time in seconds 


instantaneous distance from the satellite to the point. of 


observation, i-e., range 


anf 

angular gyromagnetic resonance frequency, 2m fpj 
angular plasma frequency, 24f, 

number density of electrons in electrons per meter} 


permittivity of free space 


ii 


e¢ pop € 


1 7 


dh 


- propagation constant (appendix only) 


polarization rotation due to the Faraday effect in radians 
length of the mean propagation path in meters. 

angle subtended at the earth's center by r 

electric flux density 


angle between the magnetic field vector and the direction 
of propagation 


element of actual propagation path length 
element of optical path length 
element of height normal to the earth's stirface 


angle between the direction of propagation and the zenith 
at the satellite 


7 -2 ee 
constant equal to 2.97 X 10 ; rationalized MKS units 
constant equal to 2.39, rationalized MKS units 


zenith angle of the propagation path at the earth 
observation point 


H cos > sec for a flat earth, H cos $ sin(i- §)/(sin i - sin §) 


for a spherical earth 

value of M at the satellite 

velocity of light in a vacuum 

critical frequency at hm 

height of the maximum of the F, layer = hmax 

mean value of M (see Appendix V) 

time rate of change of the Faraday polarization rotation 


iv 


N, = 


2 = a 
a mM mM 


eo 
iT) 


number density of electrons at the satellite height 


= magnitude of the magnetic field at the satellite height 


= scale height (usually in kilometers) 


Hmax = electron density at the maximum of the F, layer 


= angle between the magnetic field and the direction of 


propagation at the satellite 


= correction factor for neglecting higher order terms. in the 


binomial expansion for © (main body of the report only) 


ee ee 
SHinax = N dh; integrated electron density to hy, 


0 


Nyp = integrated electron density to hg computed from satellite 


positions a minute apart assuming the dN/dh ~ 0 


NH = theoretical value of the integrated electron density to hy based 


a 
" 


“ 
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on an assumed distribution N(h) 


Na/Np = ratio of ionization above F, max to that below 


= height of the equidensity column having density Nmax and 


total content Ny - SHmax * 
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As an artificial earth satel r onized 
atmosphere of the earth, the may be altered 
by its interaction with the ionosphere. The magnitude of this effect 
is related to the density of the electrons at satellite altitudes. It is, 
ectron content in the i 


ion as well as to form a 86 


omplete iono- 


sphere. 


Many methods have been used to measure the number density of 
electrons. lying above the max 
refraction,’ Faraday pol nm rotation of luna 
r gnals;* polarization rotation of artificial eart 
satellite doppler measu: ;;? radio absorpti surements,” 
yt dispersive doppler measurements , : 
: ce Administration, ’? 


um of the F, layer: radio star 


flected radio me 
e radio sig rales 
~ 0! 


lz. 


direct measurements, by National Aeronautics and Spa 


and various experiments using rockets and satellites. 4915 A rather 


plete historical record of the early work is given by Evans," who 
¢ Observation methods. All the above 


ionosphere; but the accuracy of the results is often ambiguous. This 


ty arises chiefly from two factors: (1) the ionosphere, 4 virtual 
noving particles, must usually be treated as a “quiet™ spherically 
stratified region in order to interpret results, and 


mations and assumptions must be made in analyzing 


the results (e.g-» 


approximating the ray path of a radio signal in the ionosphere). These 
nproved, 


measurement techniques, however, are continually being 
and this report is primarily concerned with some new ingi 
the method of analyzing the rate of polari 

signal in order to gain more accurate el 
upper ionosphere - 


tion rotation oO 


ionosphere by Brown’ et al.in lunar ref This 
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group set down the first order equation to determine the electron content. 
Bowhill, 17 goon after the advent of artificial satellites, reported that 
$y Constant for a 


ols* set down 


the rate of rotation is, under certain assumptio 


aaa divection of satellite motl n- meuer: and se 1é 
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ra “of pan eisation rotation fenren more e ge 
and at two ae eens Larasge 2 ‘This work was 3 
that of Arendt.” 
= have. = 
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aday pense 

! cepoded ba. ame an : srriott, 1) Little 

' and Yeh and Swenson.* These -* e a ements | have 

soak data on bid total lectron. 0 


This Laboratory has worked on the meagurement | of iOF dee apes 
electron densities since the advent of earth sat F ily 
hrow age = were ot} Hame - 26 ‘The radtio sect 
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tails of the derivation of Eq. (1) are given in Appendices 1 T ate 
II. Several iit assumptions are implicit in Eq. (!): the gy 


frequency, fu» 2 and the critical frequency, fo, are assumed soni 
smaller th nthe 6 signal frequency, th quasi- longitudinal condition 


smaller than t 
of propagation must exist (see Appendix II), and the ionospher 
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assumed to be spherically stratified. If the further assumption is 
made that the frequency f is high enough that the higher order terms 
in Eq: (1) can be neglected, then the polarization rotation is given 

( aS 

{ 
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NH cos $ ds: 


An alternative form of Eq. (2) which is useful instanteously is 


in the : 
of aed 


r density 
de, hy it is con- 


(6) 


16 3 
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in hin aa Ne ix auld, as. well as sie electron density distribution assumed in 


Since the integrated electron density is wsually of interest, 
Eq: (6) can be written as 
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(7) Nyy 
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2. The Faraday-Rate of Rotation 


The use 
limited use | because ore { 
at a single 


alter rnate em 


a ion 


= .3 Ki | ¢ oN 
(8) 1 6s ne (H eos $)ds *f H cos > x ds 
+ eT 


ds | 
8 at | 


+N gig cos > 


Another form, assuming an optical Fay path, is the derivative of 
Eq. (6), and is 


38 . Ki -_ py Oe 
(9) ras “7 dh # M Ne “ae a 


vi we Fé the ageatdapes » By sep to ane values. at the satellite. Inte- 
8 |: (10) by parts puts this equation in 


a “more night ions’ vd om 


+ + M, N eos § = ae 


ay ££2+4 | na +m, | # 


h 
° \ at 


In the derivation of Eq. (10) and Eq. (11), the fact was \ used t 
when h = ©. It is obvious that Eq. ( 11) could have been ol} 
Eq: (8) by differentiating and adding and subtracting the | 


that N= 0 
ed from 
roper terme. 


Equation (11) is a desirable form of the first order Faraday rate 
of polarization rotation since the first thi ee te: rmé can be coneide red 
to yield approximate results with the r ) 
error. The first three terms are mu 
those of Eq: (8) and Eq. (9) because ne 
under the integral in the first two terms, 
calculate and can be found more 
distribution need not be assumed. 
not necessarily small. 


¢ ip one. penrnnl re na in, 
and Mg_ ig much e 


3. 


The first two terms of Eq. (I) may be differentiated to yield a 
second order expression for the time rate of change of the polarization 


plane. Performing this operation, d6/dt is given by 


fo leg Ky ee aS | 
(12) — =“ +\ N5 — ; 1 cos $)ds + \ H cos 2 & dg + NgH, cos: oot. 


0 “o | 
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of thé last three terms. On this basis, one can cal 
slander eel of rae ae and eg aie terms of E 


ution of the 


|Fourth term} = 2.02 x 10"? Nob|First term] 


where No is the maximum electron density of the ionosphere. Assuming 
an example for noon, No = 10 electrons / cubic meter, 


(2.02 x 107" )(20"2 )(.088)| First term| 


At 


|Fourth term] 


= .0176 | First term| 


iginates 75 km above the Fz; layer maximum 
be + a 2a). 


¢€. PHYSICAL PARAMETERS AND CORRECTIONS 


TO THE THEORY 


As discussed by Mitra, ** ionized layers result from the fact 
that absorption of solar radiation penetrates an atmosphere in which 


age — increases with hase ‘ _ Thus the meets: de crea. ees 


1 ae. 
dosapted hypotheetey and y his renting exp 


= im (| } 
(13) N= No exp = {3 -Z- eZ} 
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where 
| h = height from the earth's surface 
ne hp, = height of maximum ionization 
¢ 


No = maximum number density of electrons at hyp, 


H® = kT,/mg = scale hieght 
k = Boltzmann's constant 
Ta = absolute temperature, 
m = mass of electron, and 
g = acceleration of gravity. 


eters re catch. 


Several numbers for the scale height have been found experimentally 
by assuming Chapman distributions; Yeh and Swenson® have reported 
H°=98 Km, Wright 26 has found 100 km; and Garriott® has found H® = 
65 km yields good agreement with experiment just above the maximum 
Le of the F, layer. The scale height, however; is not, in general, a con= 


Lo oe stant with height, but according to recent studies,?7 »2* the ionosphere 
is approximately isothermal above the maximum of the F2 layer. This. 
if means the scale height variation with altitude is inversely proportional 
a to the acceleration of gravity, g, in the upper ionosphere. The scale 
height variation with time, however, is proportional to the asymptotic 
[ temperature, ergo, the intensity of the solar activity. 


In previous measurements at this Laboratory,?? it was found that 
[ the Chapman distribution and a scale height increasing linearly with 
height from the F, layer maximum at a rate of 5 km per 200 km gave 
good agreement with measured results of integrated electron density. 
This is the electron distribution which is assumed in this report where 
such assumptions must be made. A typical noontime variation of the 


electron density using this. mode! is shown in Fig. 3a. 
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Fig. 3a; Chapman distribution. 
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The variation with time of Nyy, as previously stated, depends. on 
the solar activity having certain regular periods of variation: diurnal, 
seasonal, and sunspot-cyclical. A typical daily variation is shown in 
Fig. 3b which shows a two-month average on the integrated electron 
density to infinity, as calculated by Wright and Fine? 
soundings and a 100 km scale height Chapman dist 
seasonal and solar-cyclical variations may be found 


in Mitra .>? 


In addition to the time and height variations of electron density, 
horizontal ee exist in the ionosphere. These gradients often 
causic erimental reaulte* id ! 
veldedty santos to the earth's surface ig so high that t! 


quency may change along the line of tapes to the smngeee 
is a need — . pecey thet e 


discussed, it is in order to discuss some sical q 
in Eq. (11). The factor Mg, for the flat earth approximation is 
Mg 
where, to reiterate, Hg is the magnituc 
satellite, $g is the angle between the mag 
erection of PCE REIN | at the a ite, ane 


WW 


Hg cos $g sec € 


maousde field is pata to he patioks rs a 
(see Appendix IV), and since the geometry will ae ea constant for 
a fixed observation point (see Appendix III), contours for Mg can be 
drawn: Figure 22 shows the value of M, for three different satellite 
heights as a function azimuth and vertical angle from Columbus, Ohio. 


On inspection of Fig.22it ia obvious that the time rate of change 
of Mg as a satellite passes over the observation station is, for a fixed 


velocity, a function of the direction of travel and the nearness of approach. 
An infinite variety of Mg-versus-time curves is, of course, possible. 
Three typical Explorer VII passes are shown in Fig. 4 and the geometries 


are shown in Fig. 5. 


The rate of change of the path length dse/dt, for anes report is. a6- 
sumed to be the change in the optical path pai or the time rate of 
change of the actual range. This assum ahah if refraction and 
1 geometries 
ies a. 
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Equation (14) will be examined later in the light of experimental results 


to check its validity. 
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2. Corrections and Discussion. of the Theory 


Equation (2) was derived associa that the higher order terms: 
Ae could ad ee Isabel however, has outlined a 

gh-£ ea agrats assumption while 
A deed Laat 
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, fon eoa > da. 
Py) 


\- NH eos > dat = BC eo H cos $ ds * +++ 
~“0) 


iT | 


H = 75 + .025(h-h,,,), the assumed scale height, 


= 80.6 No é : 


red 
(oJ 

7 
iy 


iY 


No = the electron density at hy, and 
H = magnetic field as given in Appendix IV. 


Here, Xo, is.a convenient parameter which adapts the calculation of 

¥ toarbitrary frequencies. The calculated correction value for ¥ 
versus height was computed “et the first five terms of the binomial 
expansion for ®. The results for several values of Xo are shown in 
Fig. 7. Using Fig. 7 and ionoso onde data, ** it is possible to find the 
correction for neglecting higher order terms, and a more correct form 
of Eq: (2) becomes: 


\ NH cos $ ds 
No 


where ¥ is correct only for vertical propagation paths. 
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Chosen as discussed previously: a Chapman distribution with a scale 
height of 75 km below hy, and a scale-height gradient of +.025 above 
the maximum. It is seen in Fig. 8 that for typical Explorer VII heights 
(600 to 1100 km) the error term is much too large to neglect. It is 
necessary » therefore; at these heights to choose an electron distri- 
bution in order to evaluate 8M/dt. This procedure introduces an 
inherent error since the assumed N(h) will probably not have the same 
form as the actual N(h). Thus, the final form of the first-order 
equation relating the integrated electron density and the time rate of 
change. of the polarization plane at the point of nearest approach and 
at low zenith angles is 


alg. Pou Pa ia 7y¥ 
oN dh -=—————— 


where T iis found experimentally and is defined in terms of the rate of 
change of polarization rotation as 


3. The Two-Point High Pass Approximation 


Equation (19) will yield a number for Nyy which is subject to 
errors due to small irregularities in the ionosphere: This is an ad- 
vantage if these irregularities are of interest, but the instantaneous 
value for Njy may be sharply affected. Therefore, it is often advan- 
tagéous to average the data over a period of a minute. This is done 
as follows: 
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Using the values of T at near approach for 28 overhead passes of 
1959 Iota I, the integrated electron density, Ny, was found by Eq. (19). 
The results are shown in Fig. 11 along with the value of SHmax at 
Fort Monmouth, New Jersey. SHmax is the integrated € ) 
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Since we can only measure the magnitude 
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» Ng) Mg, and cos § are always positive, the 
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The result of studying this |\d@/dt | null position is that it yields order= 
of-magnitude information on the ratio of Ng to Ny. In Fig. 16, three 
passes illustrating three different cases are shown. It is seen that if 

is/Ny > 10°3 then the null position occurs essentially at the clos 
At of approach (CPA). As this ratio becomes less, however, the 
i iderably before the CPA, or after; depending 
h the {/ dit us, under the assumptions involved, Ny/Ns. 
ats be found from the |.48/ dt | null position, and if Nyy can be determined, 
then Ng is given directly, again, with many assumptions involved. 
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E. CONCLUSIONS AND RECOMMENDATIONS 


Values for the integrated electron densit} 
may be found from single-frequency satellite s 
of these bctaesl is. —. over our heii E 
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la. The accuracy 
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average values of SHmax as a function of tin 


late afternoon build-up in Ny as shown by Ross.” 


Experimental values of the electron content above the F2 layer 
maximum to that below show that large ratios. (N,/Np = 6 to 11) can 
C r at any time of day and not just during the evening. The average 


diurnal values, however; are slightly higher in the evening. The average 


values obtained here are generally higher than those obtained by other 
investigators using similar or other methods. 


No correlation of the experimental values of the integrated electron 


rbances 


densities to the satellite could be shown with magnetic digi 


Values of integrated ele — n density; NHP, obtained from a 
longer observation period of one minute as Opposed to six seconds, 
€ res ans ts. In fact, errors. in the satellite 
position as a function of time often showed up in these values as. 
extremely large values of Nyy. The need for more accurate time 


cor rections to satellite predictions were made obvious by this calculati 


tended to give more errati 
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Values for the electron density, Ng, at the height of the satellite, 
can bé found approximately by studying the point in time where the 
magnitude of the time rate of change of the Faraday rotation, |.d0/dt|, 
goes to zero. Ratios of Ng/ Ny > 10° yield a null position in ldevae 
at near approach. Lower values of this ratio retard or advance the 
position of the null depending on the geometry of the satellite pass. 
The accuracy of finding Ns depends upon the accuracy of finding NH 
and accounting for variations in the ray path from the satellite to the 
ground as a function of the satellites position. 


This report constitutes an improved look at the analysis of the 
Faraday effect as. applied to satellite emissions. Yet, the complete 
analysis of these satellite signals over long periods of observation is 
so complicated that this report is intended only as a preliminary 
examination of the overall problem of obtaining electron density 
information in the ionosphere from the satellite signals. Specific 
improvements that must be made concern (1) a more accurate dete r- 
mination of the satellite location with time, perhaps by using inter- 
ferometric techniques, (2) a better physical picture of the refracted 
ray peths normally encounte red as a function of ionospheric conditions 
— the B eerreeery of ‘the sere tte Labels ick a po leesiaa nee Lane hes 
represe yout tans iene Speman S jacarises eae diffex rent Spice 
frequencies since these introduce new known parameters into the 
analysis and obviate a knowledge of the satellite antenna position. 


The rate of Faraday polarization rotation of satellite signals 
provides a sensitive tool for studying the ionosphere, especially if the 
f requencies of dling OR ite Sasiot gene the me KienNS Seniery ee 
of bs aed dat re. 
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APPENDIX I - DERIVATION OF THE PROPAGATION 
CONSTANT FOR AN IONIZED MEDIUM 


This Appendix follows the previous work of Hame and Stuart.“ 
In the derivation, the medium is considered to be collision-free and 
spherically stratified with respect to the earth. 

The equation of motion of a charged particle in a magnetic field 
and an electromagnetic field is 


en by Newton's second law 


(23) Fee(E+ VX B)= mo 


Fig. 1g. Coordinat res: sy ster mM « 


direction of pr opagation : r of the ele ectromagnetic wave chosen to lie in 


the x-z plane. 


O in’ Vv oe 2 2 
22 On suming that | E and V are of the form E = E,e ef and 
V = Voes, Eq. (23) becomes sacs 
(24) 
which yields for the velocity components: 


nt 


(25) 4 


The above Eq. (25) shows that Vy, and Vy become infi 
denominator goes to zero. This angular frequency 


(26) wy = <5 


is. the sera cet frequency and this resonance Occurs, for typucal 
h's jie rie baieade at +n = 00 eps 


s. for } 
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Maxwell's first curl equation in a regi 
electrons 


On COR ntaining moving 


may be written as 


whe re 
fi e 
| Sn JEag 
(30) jk | j€h2 €yy 
I 0 (0) ry 
Equation (30) defines the tensor pe rmittivity or dielectric constant 
dyadic which resuits pa difference in direction of E and D. This 


41 


tensor has been (presented wd bpd poms oa Whit tmer. 


where f, denotes the plasma frequency . 
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Maxwell's second curl e quation for the , g given time dependence isi 


On substituting Eq- (32) into Eq. (29) the following equation is ob 
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where E, for our purposes, has the form 
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From Fig. 18, the propagation constant Y has < 
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(35) Y= ¥yX FF ¥Z2 5 Fy, Xand Zz unit vectors. 
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always ak and in this case ered ‘coeff ient dete rove is ZETO. 


Using the fact that the determinant is zero, Y can be solved for 
i lements. The result is 


: = ( 2 a2 

+2:€ a €33, + (em = E12 = =€qy € 3 
{ 2 V2 

+4655 €29 cos? : : 


2 
head = ety =€qy €33 


2(€14 - € 33, ) sin* > + 2€ 55 


a = 3 | ¢ 
(40) ¥' il el] 
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which is the expression for the propagation constant of electromagnetic 
waves inan ionized media suiting the assumptions involved in the 


derivation. 
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sless, the propagation 
constant may be replaced by the phase constant: 
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the polarization plane per unit distance: 
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represents the angular increment of fatok 
polarization over a distance dr. 
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